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By use of the spectral density approach the influence of 
the lattice structure on the possibility of ferromagnetism in 
the single band Hubbard model is investigated. The d = oo 
hypercubic lattice does not show magnetic phase transitions 
of second order irrespective of the strength of the Coulomb 
coupling. However, first order transitions to finite magnetic 
moments, not visible as singularities of the paramagnetic sus- 
ceptibility, may appear in the very strong coupling regime. 
In d=3 second order transitions are found but only for very 
strong couplings, where the non-locality of the electronic self- 
energy acts in favour of the spontaneous magnetic moment. 
The influence of the non-local part of the self-energy is partic- 
ularly strong for lattices with small coordination number. The 
non-bipartite fee lattice exhibits saturated ferromagnetism for 
all band occupations 1 < n < 2 while for less than half filled 
bands (0 < n < 1) the system remains in any case param- 
agnetic, and that for d = 3 as well as d = oo. The Curie 
temperature runs through a maximum at about n = 1.4 and 
vanishes for n — » 1 and n — » 2. 



I. INTRODUCTION 

Originally the Hubbard modelEH! was introduced to 
describe the bandmagnetism of itinerant electrons, at 
least qualitatively. The model Hamiltonian 



W=£(Ty-M)4 



i,a 



(1) 



incorporates, in the simplest way, the interplay of kinetic 
energy, Coulomb interaction, band structure and Pauli 
principle for studying magnetic and electronic properties 
of interacting electrons in a single non-degenerate narrow 
energy band. is the intersite hopping integral, usually 
taken between nearest neighbors, while U represents the 
Coulomb repulsion. c$ CT (c ia ) stands for the annihilation 
(creation) operator of an electron with spin a =|, j at 
lattice site R^. By [i we denote the chemical potential. 
rii a = <^ ia Ci a is the occupation number operator. The 
model-Hamiltonian defines a non-trivial many-body 
problem that could be solved up to now only for some 
limiting casescTTj. 

Although it was soon realized, that the Hubbard model 
is rather a generic model for antiferromagnetism, the 
question whether or not the Hubbard model possesses 



ferromagnetic solutions continues to be the matter of 
controversial discussions. In this work we focus exclu- 
sively on ferromagnetic solutions of the single band Hub- 
bard model. What concerns ferromagnetism of real sub- 
stances, as Fe, Co, Ni, the direct Heisenberg exchange, 
Hund's rule coupling ^nd band-degeneracy, of course, 
play an important roleta, the investigation of which is be- 
yond the scope of this paper. 

The Merinin- Wagner theorema applied to the Hub- 
bard modell3 excludes ferromagnetism for one- and two- 
dimensional lattices. For the three dimensional system 
with infinitely strong Coulomb repulsion (U — ► oo) ex- 
act results have been derived in a pioneering paper by 
NagaokaO. For very special band fillings n± = (N±l) /N 
(N: number of lattice sites) ferromagnetic ground states 
are possible; in the case of bipartite (sc and bec) lattices 
for n + as well as for n_, in the non-bipartite lattice (fee) 
only for n + . This indicates that the possibility of ferro- 
magnetism crucially depends on the lattice structure. For 
the hynercubic (he) lattice at infinite dimensions Fazekas 
et. al.U showed that saturated ferromagnetism is ex- 
cluded. However, in a recent study UlmkeEZl applied a 
finite temperature Quantum Monte Carlo calculation to 
the Hubbard model on a fee lattice at infinite spatial di- 
mensions. A highly polarized ferromagnetic ground state 
is found in a rather wide range of band occupation. It 
is to be supposed that the greatest chance of ferromag- 
netism to appear in the._sinp;le band Hubbard model is 
for non-bipartite latticesESIiJ. 

In this paper we want to investigate how sensible band 
ferromagnetism in the Hubbard model depends on the 
lattice structure, and therewith on the density of states 
of the non-interacting electrons (Bloch density of states: 
BDOS). We apply a "spectral density approach" (SDA) 
the reliability of which, at least for the strong coupling 
limit, has been demonstrated in previous papersEjliZl. 
The same approach has been applied to a somewhat gen- 



ubbard model) 
. The SDA ba- 



eralized model of magnetism (multib 
for the band ferromagnets Fe, Co, Nil 
sically consists in a two-pole ansatz for the single-electron 
spectral density, justifiecL-by the rigorous moment anal- 
ysis of Harris and LangaHJ. The free parameter of the 
ansatz are self-consistently fixed by equating exactly cal- 
culated spectral moments. The SDA tums-aut to be es- 
sentially equivalent to thc_Ecith mcthocciEJ and to the 
Mori projector formalismeJnJ. 

For several types of BDOS we have derived the static 
paramagnetic susceptibility x as a function of the band 



1 



occupation (0 < n < 2) and the temperature T. From 
the singularities we read off the instabilities with respect 
to ferromagnetic order. The results are confirmed by a 
direct calculation of the spontaneous magnetic moment. 
However, the x _1 -zeros are due to magnetic phase tran- 
sitions of second order only, while additional first order 
transitions to a finite spontaneous moment may occur. 

The SDA statements on the possibility of ferromag- 
netism in the Hubbard model can be compared with re- 
cent results from the d = c»-techniquec3iij. As to the 
lattice structure dependence the SDA solutions predict 
qualitatively the same trend. It is found, e.g., that for the 
d = oo-hypercubic lattice ferromagnetism is excluded, 
being, on the other hand, possible for an d = oo-fcc lat- 
tice. 

For finite spatial dimensions (d = 3) we find that the 
non- locality of the quasiparticle self-energy, vanishing for 
d — > oo in the SDA, too, has a stabilizing influence on 
the magnetic order. As an example, the d = oo hy- 
percubic lattice is non-magnetic, but the d = 3 sc lat- 
tice becomes ferromagnetic for U/W > U c /W = 4 and 
0.34 < n < 1.66 if the full k-dependent self-energy is 
taken into account. Neglecting the k-dependent part of 
the self-energy drives U c /W to a higher value (= 14) with 
n c only slightly changed. 



II. SPECTRAL DENSITY APPROACH 

The method is based on a physically motivated ansatz 
for the single-electron spectral density, which is defined 
by: 



Si ja (E) = J dte-t Et S lja {t) 

— OO 



(2) 



(3) 



[.. , ..] + is the anticommutator, and (. . .) means thermo- 
dynamic averaging. The operators in (j^) are thought 
as time-dependent Heisenberg operators. The mentioned 
ansatz contains some free parameters which can be fitted 

(n) 

by equating a set of spectral moments M-jJ 

+oo 

?i = 0,1,2, 



M 



(n) _ 



ij a 



JdEE n S lja (E): 

— oo 

The moments are calculated exactly via 
. . [c irT ,H]-, 



M (n > 
13° 



.n 



(4) 



(5) 



n-fold commutator 

In ref.0 a two-pole ansatz for the spectral density of the 
strongly coupled Hubbard model is justified as a reason- 
able starting point requiring the fitting of the first four 



spectral moments. The evaluation leads to an electronic 
self-energy of the following structure: 



^ A (E) = U{n _ a ) 



E + n - B. a - F k _ 



E + ix - B_ a - F k _ a - U (1 - (n_ a » ' 



(G) 



The decisive terms arc B- a and F^ a . For rB- a = 
Fk-CT — @ reproduces the Hubbard-I solutiontl, which 
is incapable to describe ferromagnetism. -B_ CT as well as 
Fk-cr consists of "higher" correlation functions, possibly 
creating a spin-dependent shift of the Hubbard bands 
and therewith resulting in a finite spontaneous magnetic 
moment. The term B_ a , 

, i& 

(n-a)(l - {n-a))B- a = — ^Ti^cl^Cj-vin^ + n ja - 1)}, 

id 

(7) 

can rigorousl y h e expressed by the single-electron spec- 
tral densityEJO. The other term, 



(n_ ff )(l-(n_ CT ))F k _ ff 



3 -ik-(R,-Rj 



N 



{(ni-rUj-o) - n 2 _ a ) 



incorporates a density-density, a double hopping-, and a 
spinflip-correlation. For the evaluation of these expecta- 
tion values properly modified spectral densities, 



S ( i fJ(t) = ±([A i (t),cl(0)} + ), 



(9) 



are introduced, where the operator A, is chosen in such a 
way that the respective expectation value can be derived 
directly via the spectral theorem from S^J(E). It can 
be justified by inspection of the Lehmann representation 
of these functions that in the wave- vector representation 
they must represent two-pole functions as the original 
single-electron spectral density The various func- 

tions (||) and (^|) differ only by the spectral weights of the 
poles. For fixing S^J (E) one therefore needs only two re- 
spective moments. Doing so one finally arrives at a closed 
system of equations that can be solved self-consistently 
for all quantities of interest. For further details the reader 
is referred tcO. 

We study the magnetic phase transition by use of the 
static susceptibility 



x(T,H,n) 



1 



(dM 
Mo \~dH 



(10) 



fiQ is the vacuum permeability and H a homogeneous 
static magnetic field. M denotes the magnetization 



2 



N 

M=- MB ((n T )-(n;». 



(11) 



We calculate \ for the paramagnetic system in the zero- 
field limit. The poles of the susceptibility indicate the 
instabilities of the system towards ferromagnetic order- 
ing. To calculate the static susceptibility x we have to 
include a Zeeman term in the Hamiltonian (Q). When 
differentiating M with respect to the external field one 
has to bear in mind that all expectation values (...) are 
field-dependent. 



III. RESULTS 

We have investigated the tendency to ferromagnetism 
in the single band Hubbard model for several lattice 
structures and different spatial dimensions. The zeros 
of the inverse static paramagnetic susceptibility indicate 
the magnetic phase transitions of second order. Fig. [l] 
shows the results for a sc lattice. 
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FIG. 1. Inverse paramagnetic static susceptibility \ f° r 
the sc lattice (d = 3) as a function of the band occupation 
n for various values of the Coulomb interaction U. (a) Sys- 
tem with the full k-dependent self-energy, (b) System with 
a local self-energy (Fw-a = artificially, see eq. (^|)). The 
inset in (a) shows the magnetization curves m(n) of the two 
ferromagnetic solutions starting at the two zeros of y _1 for 
U/W = 5. (T = OK). 

Sometimes it is arguecfi that for this bipartite lattice 
ferromagnetism is hardly to be expected. On the other 
hand, Shastry et. al.EZl report for the d — 3 model even 
saturated ferromagnetism for n > 0.68. That is very-dose 
to our findings, ft is a special feature of the SDAEjIIZI, 
possibly even for the Hubbard model itself, that there 
are two ferromagnetic solutions, i. e. two zeros of x^ 1 ■ 
The first solution sets in at — 0.34 (U — > oo), where 
the actual value only slightly depends on U. This so- 
lution runs into saturation for n > 0.68 (U — > oo), in 



exact agreement with the results of ref.0. The second 
solution appears for higher band occupations, but does 
never reach saturation. This less magnetized solution is 
always less stable and can be disregarded in the following 
discussion. In spite of the appearance of ferromagnetic 
solutions the sc structure seems not to be very conve- 
nient for a spontaneous order. A rather strong Coulomb 
coupling U/W > 4 (W: Bloch bandwidth) is needed. In 
this context it is interesting to look at the influence of the 
non-locality of the electronic self-energy (^) . The wave- 
vector dependence comes into play due to the higher cor- 
relation functions in Fk_ CT (||)- Since in the limit of infi- 
nite spatial dimcnskms the electronic self-energy is wave- 
vector independent^, the role of the k-dependent terms 
in Y^ A (E) might be underestimated for d = 3. In part 
(b) of Fig. [l] we have plotted x _1 a $ function of the elec- 
tron density for the case that the self-energy @ has been 
made local simply by setting Fk- CT = 0. The critical 
couplings U c /W for ferromagnetic order rises enormously 
from 4 for the full k-dependent self-energy (Fig. 0(a)) to 
14 for the local self-energy (Fig. 0(b)). Its influence on 
magnetic stability, however, is not for all cubic lattices 
of the same importance. It can be recognized that with 
increasing number of nearest neighbours the importance 
of the k-dependent part of the self-energy rapidly de- 
creases. For theJjcc lattice [Z\ = 12) it appears already 
non-remarkabfyli3. 

In infinite spatial dimensions, for which the higher cor- 
relation Fk-o- disappears, we do not find zeros in x^ 1 f° r 
the hypercubic lattice (inset of Fig. ^) , and that irrespec- 
tive of the strength of the Coulomb coupling U. Phase 
transitions of second order are therefore excluded. 
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FIG. 2. Magnetization m for the hypercubic lattice in 
d — <x as a function of the band occupation n for various 
values of the Coulomb interaction U. The BDOS is given 
by: po{E) = exp(-0.5(£/t*) 2 )/(t* V2tt). Inset: Inverse para- 
magnetic static susceptibility x _1 as a function of the band 
occupation n. (T — OK). 
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The same fact has been reported by JarrelEa who elab- 
orated a self-consistent Monte Carlo procedure, which 
is claimed to be essentially exact for the d = oo Hub- 
bard model. Although there is no indication for ferro- 
magnetism in the static susceptibility, the direct evalu- 
ation of the SDA for the spin-dependent, averaged par- 
ticle numbers (n-j-), (n{) reveals ferromagnetic solutions 
(Fig. ||). They belong to first order transitions being 
therefore not visible as x _1 - zeros - These solutions ap- 
pear for rather strong Coulomb couplings and for band 
occupations n > 0.57. For a given U the m(n) curve 
forms a closed bubble steadily contracting itself for in- 
creasing temperature. For a given band occupation n this 
means a first order transition at a certain critical tem- 
perature. Ferromagnetic saturation is, strictly speaking, 
never reached, not even for T = OK (Fig. ||). This is 
because of the infinite tails of the Gaussian BDOS. The 
strict physical meaning of this kind of ferromagnetism in 
the Hubbard model requires further investigation. 

It turns out that the three dimensional bec lattice with 
its divergence at the center of the BDOS is much more 
convenient for ferromagnetism than the sc lattice. Al- 
ready for U/W ~ 1 two ferromagnetic solutions appear, 
the stable one starts for U — > oo at n c = 0.52 reaching 
saturation (m=n) for n > 0.68 (Fig. 0). 
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FIG. 3. (a) Inverse paramagnetic static susceptibility x 
for the bec lattice as a function of the band occupation n for 
various values of the Coulomb interaction U (parameters as 
m Fig. 0(b)). The thin dotted line corresponds to the system 
with the local self-energy (Fk-a = artificially, see eq. (|^)) 
for U /W = 5. The inset shows the magnetization curves m(n) 
together with x _1 f° r U/W = 5. (b) the same as in (a) for 
the fee lattice. (T = OK). 

The latter result is again in exact agreement with the 
findings in ref£3, where the stability of the Nagaoka 
statetl (U/W — > oo) with respect to an electron spin flip 
is investigated. As already mentioned the k-dependencc 
of the electronic self-energy has not such dramatic con- 



sequences as for the sc lattice (Fig. [j]). 

The non-bipartite fee lattice does not allow ferromag- 
netism for less than half-filled energy bands. However, 
for the more than half-filled band we find saturated fer- 
romagnetism (to = 2 — n: "Nagaoka state") for all elec- 
tron densities 1.0 < n < 2.0 indicated by a zero of x 1 
at n — 2.0 (Fig. ^). This_agrees with the investigation of 
Mullcr-Hartmann et. alO who employed the single, spin 
flip Gutzwiller variational states of Shastry et. al£3 and 
corrected their value of n c = 1.62 ir£3 to n c — 2.0. For 
electron densities 1.0 < n < 1.2 the SDA yields a second 
ferromagnetic solution which is less stable than the fully 
polarized one. 




! FIG. 4. fee lattice 

(BUOS: po(E) = exp(-0.5(l + V2E))/ y/-K{l + V2E) from 
ref.li3) (a) Inverse paramagnetic static susceptibility x~ 1 ( n ) 
as a function of the band occupation n for various temper- 
atures T (solid: T = OK, long dashed: T = 500 K, dotted: 
T = 1000 K, dot dashed: T = 1500 K, dashed: T = 2000 K, 
dot dot dashed: T = 2500 K.) The inset shows the Curie 
temperature Tc as a function of the n. The solid line is the 
result of the SDA, the long dashed line is taken from ref.E3. 
(U — 4eV). (b) x _1 ( n ) for various values of the Coulomb in- 
teraction U. The inset shows the magnetization curves m(n) 
together with x _1 f° r U = 4eV. (long dashed: U = 0.4 eV, 
dashed: U = leV, dotted: U = 4eV, solid: U = 100 eV; 

t = ok). 

The results are qualitatively similar for the Hubbard 
model on the infinite dimensional fee -Lattice (Fig. [!]). 
Here, we used the BDOS given in ref.LL3. The x _1 (n) 
curves exhibit two zeros due to two ferromagnetic solu- 
tions as in the d — 3 case (Fig. ||). The fully polarized 
state (m = 2 — n) is stable for all band occupations 1.0 < 
n < 2.0 with a critical coup ling U c (n = 2.0) /W = + . 
The same is reported in ref.E3li-X The asymmetric BDOS 
of the non-bipartite fee lattice is obviously very conve- 
nient for the band ferromagnetism in the more than half 
filled single band Hubbard model. 
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For finite temperature the \ 1_ zero shifts away from 
n = 2.0 to lower values indicating a Curie-temperature 
of T c = 0+ for n = 2 (Fig. |(a)). The T c ( n) curve runs 
through a maximum at n R* 1.4 and goes to zero for the 
half-filled (n — » 1.0) and the fully occupied (n — > 2.0) 
energy band (see inset in Fig. |](a)). Thes^. results are 
qualitatively ccpsistent with that of UlmkeO. However, 
contrary to ref.113, our Tc(n) curve persists until n = 1. 
This reflects the fact, that close to half filling (1.0 < 
n < 1.15) the magnetization curves in the SDA exhibit 
first order phase transitions as a function of temperature. 
Magnetic solution of this kind where not considered in 
refJ13. Quantitatively our X^-values are on an average 
higher by a factor 4, but nevertheless in a reasonable 
order of magnitude. 



IV. CONCLUSIONS 

By use of a spectral density approach we have demon- 
strated the striking influence of the lattice structure on 
the possibility of ferromagnetism in the Hubbard model. 
The zeros of the inverse static paramagnetic susceptibil- 
ity indicate the onset of ferromagnetic order in depen- 
dence of typical variables as the band occupation n and 
temperature T. In the he lattice no "normal" ferromag- 
netism appears at infinite dimensions, while for d — 3 (sc 
lattice) magnetic solutions are found if U/W exceeds a 
rather high critical value. For small lattice coordination 
number (e.g. sc lattice) the non-locality of the electronic 
self-energy favours ferromagnetism in a remarkable man- 
ner. In the non-bipartite fee lattice no magnetism exists 
for less than half-filled bands, while for 1.0 < n < 2.0 
saturated ferromagnetism is found at d = 3 as well as 
d = oo. 
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